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Abstract 

Nuclear pairing correlations are known to play an important role 
in various single-particle and collective aspects of nuclear structure. 
After the first idea by A. Bohr, B. Mottelson and D. Pines on sim- 
ilarity of nuclear pairing to electron superconductivity, S.T. Belyaev 
gave a thorough analysis of the manifestations of pairing in complex 
nuclei. The current revival of interest in nuclear pairing is connected 
to the shift of modern nuclear physics towards nuclei far from sta- 
bility; many loosely bound nuclei are particle-stable only due to the 
pairing. The theoretical methods borrowed from macroscopic super- 
conductivity turn out to be insufficient for finite systems as nuclei, in 
particular for the cases of weak pairing and proximity of continuum 
states. We suggest a simple numerical procedure of exact solution of 
the nuclear pairing problem and discuss the physical features of this 
complete solution. We show also how the continuum states can be 
naturally included in the consideration bridging the gap between the 
structure and reactions. The path from coherent pairing to chaos and 
thermalization and perspectives of new theoretical approaches based 
on the full solution of pairing are discussed. 



1 Introduction 



Nuclear pairing is one of the main and longstanding pillars of current under- 
standing of nuclear structure. Pairing provides an important contribution to 
the odd-even mass difference in the phenomenological mass formulae |1 . As 
an empirical fact the pairing was put in the foundation of the shell model 
by Mayer and Jensen [2 in order to be able to predict ground state spins 
and other properties of non-magic nuclei. In the shell-model framework, the 
classification of paired states is usually performed with the aid of the senior- 
ity scheme jHJ Ej , where the seniority counts a number of unpaired particles; 
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a similar scheme is used in atomic spectroscopy The Bardeen-Cooper- 
Schriffer (BCS) microscopic theory of superconductivity elucidated the 
main features of the ground state, excitation spectrum, transition probabili- 
ties and phase transition in a Fermi-system governed by the attractive pair- 
ing. Immediately after that A. Bohr, Mottelson and Pines pointed out [7] 
the similarity between the superconducting pairing correlations and observed 
pairing effects in nuclei. The thorough application of the BCS approach to 
the nuclear problem was done by Belyaev in his seminal paper jSj. It was 
quantitatively demonstrated that the pairing correlations influence nearly 
all phenomena in low-energy nuclear physics: binding energy, single-particle 
spectra, transition probabilities, collective vibrational modes, onset of defor- 
mation, rotational moment of inertia, level density and thermal properties. 

The BCS theory, as well as its advanced form the Hartree-Fock-Bogoliubov 
(HFB) method 0, is formulated in a way fully appropriate for macroscopic 
quantum systems; in fact it gives an asymptotically exact solution [TU] in the 
thermodynamic limit. For mesoscopic systems, as nuclei, atomic clusters, 
quantum dots, Fullerenes, nanotubes or small metallic grains, this approach, 
although qualitatively reflecting the main physical features, turns out to be 
insufficient. The total number iV of particles is preserved in this method 
only on average. Since we have to describe the spectroscopy and reactions 
for a specific nuclide, we need either to add special projection procedures 
[TT1 IT^l IT31 HI] that fix the exact value of N or generalize the formalism by 
approximately including the matrix elements restoring the particle number 
conservation [13 [EJ EH EE1 ■ 

Another drawback of the BCS or HFB approaches is a sharp phase transi- 
tion as a function of parameters or temperature. As pointed out by Belyaev 
jS|, in a system with a discrete single-particle spectrum, the Cooper phe- 
nomenon requires, in contrast to a macroscopic Fermi-gas, a certain mini- 
mum strength of pairing attraction. For a weaker pairing, the mean field 
approaches, as BCS or HFB, give only a trivial normal solution, while in 
reality the effects of pairing correlations still exist. The pairing correlations 
in the mean field framework also vanish immediately after the thermal phase 
transition. These predictions are incorrect for mesoscopic systems. The ex- 
act shell model calculations show that the pairing correlations do not 
disappear at the BCS transition point revealing instead a long tail of "fluc- 
tuational superconductivity" . 

The main field of interest in nuclear structure is currently shifted to the 
nuclei far from stability. As we move to loosely bound systems, the influ- 
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ence of the continuum becomes exceedingly important. Along with that, all 
attractive correlations are to be taken into account properly in order to de- 
termine the position of the drip line. The correct treatment of pairing as 
the main attractive part of the residual interaction is absolutely essential for 
such problems. Some nuclides, like the notorious n Li, are bound just due 
to the pairing correlations between the outermost neutrons, an example of 
a real Cooper pair [20j. The theory of pairing including both discrete and 
continuum single-particle levels is still in its infantry [2T] . 

Finally, there is a clear necessity to understand the interplay of pair- 
ing with other parts of the residual nuclear interactions going beyond the 
mean field approximation of the HFB method. Of course, in the lower part 
of the nuclear chart (p, sd and pf shells) there are well developed modern 
shell model methods and reliable effective interactions of the nucleons in 
the truncated single-particle space, see for example [221 123 I2H I2H] ■ With 
the possibility of incorporating additional stochastic and statistical elements 
[2EI 123 I2H| > the shell model calculations are able to describe an impressive 
amount of spectroscopic data. Unfortunately, the qualitative interpretation 
of results obtained by the large-scale shell-model diagonalization in terms 
of simple physical models is getting quite difficult as the matrix dimensions 
approach the limit of current computational strength. In addition, one needs 
to mention that the continuum problem is not solved in the standard shell 
model approach based on the discrete spectrum. Therefore the gap between 
the shell model for nuclear structure and the reaction theory is widening. 

In this situation it is alluring to first separate the pairing part of the nu- 
clear interaction and to solve the corresponding many-body problem exactly. 
As was shown in [2H] , the exact solution is numerically simple and eliminates 
all drawbacks related to the BCS approximation. At the same time, it is still 
close enough to the standard images of nuclear structure. This exact solu- 
tion can serve as a zero-order step or a background that allows one to look 
for the new approaches and approximations to the full problem, effects of 
other interactions, inclusion of the continuum, relation to the reaction cross 
sections and so on. 

In what follows we start with sketching the traditional approaches and the 
exact solution of the pairing problem. We compare the exact results with the 
BCS approximation, both for the ground and excited states, demonstrate the 
possibility of including the continuum physics and consider chaotic aspects 
of pairing, a topic practically unexplored in the literature. We complete the 
paper with the discussion of the perspectives of new approximations based 



3 



on the exact pairing solution. 



2 Approaching the solution of the pairing prob- 
lem 

2.1 Pairing Hamiltonian 

We formulate the pairing problem in the restricted single-particle space of 
fermionic orbitals assuming the Hamiltonian 



H p = ^eiojai - ^^G 12 p\p 2 . (1) 

1 1,2 

Here the subscripts 1, ... run over the complete set of orthogonal single- 
particle basis states, and we assume the Kramers double-degeneracy of time- 
conjugate orbitals |1) and |1). The pair creation, p\, and annihilation, p 1 , 
operators are defined as 

p\ = a\a\, pi = ajai, (2) 

and the double time- reversal acts as |1) = — 11). In the important case of 
spherical symmetry of the mean field that supports degenerate orbitals \jm) 
with energies ej, angular momentum j, and projection j z = m, the pairing 
Hamiltonian can be conveniently written as 

H P = J2 e A - E G n> L ) L fi ( 3 ) 

3 j j' 

where we use the operator notations, 

Nj = hj rn = Oj m aj m , (4) 

m m 

for the occupancy operators and (aj — > a,^ = (— ) J_m aj_ m ) 

Lj = \ J2H j - m *3-ma jm , L) = \ J2(-ySm*lm, (5) 

m m 

for annihilation and creation operators of pairs in a state with certain quan- 
tum numbers of total angular momentum and its projection, J = M = 0. 
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Essentially the same Hamiltonian (JIJ) can describe the situation in the de- 
formed mean field. 

The pairing interaction is defined in terms of real matrix elements of 
diagonal pair attraction, G\\ > 0, and off-diagonal ones, G12 = G21, for pair 
transfer between the orbitals 1 and 2. We limit ourselves here by the pairing 
of identical particles (isospin T = 1) in states with zero total pair angular 
momentum; the consideration of the T = proton-neutron pairing would 
require a non-zero spin of the pair. 

2.2 BCS approach 

In the BCS theory, the ground state |0)bcs of the paired system with Hamil- 
tonian Q is determined by minimizing the ground state energy with a trial 
wave function 

|0)bcs = n 1>0 («i-^pl)|o} = n 1>oUl exp (-^rfj 1°) > ( 6 ) 

where the variational parameters for each pair of time-conjugate orbitals, u\ 
and Vi, can be taken as real numbers subject to the normalization uf+vf = 1. 
In our discussion we denote \N; s . . . ) as the lowest in energy iV-particle state 
with quantum numbers s, . . . . We use the notation |s)bcs for the BCS state 
with s quasiparticles that has an uncertain particle number; for this reason TV 
is not shown, however it is assumed that the state corresponds to an average 
particle number N. The restriction 1 > means that the time-conjugate 
orbitals are not counted twice. The exponential form of the variational wave 
function shows that this state is generated as a coherent state of fermionic 
pairs; this feature can be put in a foundation of methods going beyond the 

bcs muni- 

The variational solution is given via the occupation amplitudes 

where e\ = ^e\ + A\ is quasiparticle energy and Ai is the BCS energy gap. 
The gap equation arising from the minimization of energy is 

A: = \ BCS<0| X>224|0>BCS = E T 1 A2 ' 
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The ground state as well as low-lying excited states of paired systems can be 
classified introducing quasiparticle creation operators 

a[ = U\a\ — V\a\ . (9) 

This transformation is canonical due to the correct normalization of u and 
v. The BCS vacuum |0)bcs in © can be defined as a|0)Bcs = 0. The 
Bogoliubov transformation Q mixes particle and hole states. 

One of the problems in the BCS application to small systems is particle 
number non-conservation which follows from the form of the wave function 
(JUJ). A common practice to fix the number of particles is to introduce a 
chemical potential \i by the shift of single-particle energies e — > e — \x. The 
right total particle number is restored in average through the condition 

J2ni(fi) = N. (10) 
l 

The uncertainty of the total particle number is given by 

(AiV 2 ) =W -N 2 = ^2 ni (\ - ni ) . (li) 

i 

Practically this fluctuation is of the order a/ (AiV 2 ) m 2 . Given that relative 
fluctuations go down with increasing N , the BCS solution is asymptotically 
exact in the thermodynamic limit of macroscopic systems [TO] . The con- 
tribution to the energy from the fluctuation is quadratic in the number of 
particles jH2j and thus gives an extra correction of the monopole type. Spe- 
cial methods, such as Lipkin-Nogami [TT] IT2"] techniques, were invented in 
order to suppress the particle number fluctuations. 

The particle number violation in the BCS is an example of spontaneous 
symmetry breaking with respect to phase rotations generated by the number 
operator, U(<f>) = exp(— i(f)N/2) . While the pairing Hamiltonian is invari- 
ant under this rotation, [H p , U((f))] = 0, the trial ground state |0)bcs has 
a preferred orientation, with the usual consequences of the appearance of 
Anderson-Goldstone-Nambu modes whose properties are influenced by the 
mesoscopic nature of the nuclear systems [33]. A number of projection tech- 
niques have been developed within the framework of the symmetry violation 
treatment, see [34J and references therein. 
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2.3 Recursive method with particle number conserva- 
tion 

Instead of introducing the pair condensate of the pairs with an uncertain 
particle number, another possibility was explored (331 ESI HZj, where the 
matrix elements of relevant operators explicitly keep memory of the exact 
particle number. The gap is defined now as a matrix element of the pair 
annihilation operator between the ground states |iV;0) of the neighboring 
even systems, 

Ai (iV) = l - (N - 2; 0| J2 Gi2P2\N; 0) . (12) 

2 

Similarly, the particle number dependence enters the single-particle transition 
amplitudes between adjacent even and odd systems, 

V\ (N) — (N — 1;1 \a\ \N\ 0), u x (N) = (N + 1; 1 \a\ \N; 0) . (13) 

Here one needs to consider a sequence of ground states \N; 0) with energies 
E(N) . It is assumed that the spectra of adjacent odd nuclei start with 
energies E(N±l; 1 ) of the states |iV±l; 1 ) containing one unpaired nucleon 
with quantum numbers 1. 

The exact operator equations of motion for the single-particle operators 
a\ and a{ , 

[a 1 ,H] = e 1 a 1 + - ^ G 12 a\p 2 , (14) 

2 

[a\, H] = -eia\ - G 12 ptai, (15) 

2 

can be used to construct recursive in N equations for the gap (|T2"j) and single- 
particle transition amplitudes (|13j). The approximation of no condensate 
disturbance by an extra particle jH3 EH] , 

(iV-l;l|^Gi 2 otp2|iV;0) 

2 

w {N - 1; 1 | at \N - 2; 0)(N - 2; 0| ^ G 12 p 2 \N; 0) = 2Ai (N) u £ (N - 2), 

(16) 
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leads to the recursion relation connecting adjacent even nuclei, 



v x (iV-2)| 2 = l 



A!(iV)| 2 



(17) 



[e x (iV) - e\ (iV)] 



where the A^-dependent chemical potential is introduced, 




(18) 



and quasiparticle excitation energy is defined as 




(19) 



with shifted single-particle energies 



e' l {N) = e 1 



(20) 



2 



The analogs of the number conservation equation and the gap self-consistency 
condition now read 



The pairing problem formulated in this manner allows a recursive solution 
in both directions, starting from the empty shell or from the completely filled 
shell. This solution reduces to the BCS under assumption that the gap does 
not change in the transition from N to N — 2, the same approximation of 
particle number uncertainty that lead to the BCS particle number fluctuation 
((TTJ). Based on this feature, the BCS energy can be efficiently corrected by 
the substitution A" — > A" — 1 [SZ|. Corrections to such iterative methods via 
inclusion of pair-vibration excitations in the intermediate states of Eq. (fTHj) 
with further diagonalization are also possible [TZj, as well as the treatment 
of the excitations with RPA techniques jHHl EH] ■ 

The particle conserving treatment does not resolve another problem of 
the BCS solution, namely the sudden disappearance of pairing correlations 




(21) 



where f2 is the total capacity of fermionic space, and 




(22) 
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when coupling becomes too weak. The gap equations, (jSJ) and have 
only trivial A = solutions if the pairing strength G is too small compared 
to the single-particle energy spacings. The point of this phase transition is 
roughly at the critical coupling strength G c , 



where both the pairing strength and the density of single-particle states z/p 
are taken at Fermi energy. Many nuclear systems in the shell model pic- 
ture are close to or even below the point of the BCS instability, although 
the pairing correlations do still exist |29|. As will be discussed later, near 
the phase transition, in the so-called pair-vibrational regime, the fluctuations 
drive pair scattering to an almost chaotic level leading to a sharp increase of 
the mixing between the states of the same seniority. This randomness makes 
the approximation (j!6|) or any truncation of states mixed by the pair vi- 
brations inappropriate. Various projection techniques also seriously suffer in 
the region of weak pairing. More advanced approaches, such as HFB+RPA, 
break down in the vicinity of the phase transition, though the pairing solu- 
tion can still be continued into the region beyond the critical point using the 
RPA based on the Hartree-Fock solution for a normal state. This treatment 
drastically improves the prediction for the ground state energy jJUj- The 
methods of equations of motion [321 IH] and variational techniques can be 
used to better account for the RPA ground state correlations [I2|. Being ap- 
plied to superfluid Fermi systems these methods demonstrated a considerable 
improvement |4*2*1 143j . 

3 Exact solution of the pairing problem 

Historically, a few suggestions were put forward for the exact solution of the 
pairing problem. The Richardson method, described in the series of papers 
|4*4*1 145j . provides a formally exact way for solving the pairing Hamiltonian. 
This method reduces the large-scale diagonalization of a many-body Hamil- 
tonian matrix in a truncated Hilbert space to a set of coupled equations 



for unknown parameters z\ , their number being equal to that of valence 
particle pairs. The ground state energy is then equal to E(N) = Ylfx^ z a • 



G c u F = 1, 



(23) 



(fy = 2j + 1) 




(24) 
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Recently this solution was revived and reinterpreted with the aid of 
the electrostatic analogy, similar to that used by Dyson in his theory [17] 
of random level ensembles. Unfortunately, the Richardson solution is only 
valid for the constant pairing force, Gjf = G = const. It also requires serious 
numerical efforts rapidly growing with the number of particles. Recently, 
exact solutions have been also approached with sophisticated mathematical 
tools as infinite-dimensional algebras [IB]- Such formally exact solutions have 
a certain merit from a mathematical point of view and might be useful for 
developing simple models EH] • However, they are not very promising for 
practical problems in nuclear physics. 

The natural way of solving the pairing problem is related to the direct 
Fock-space diagonalization. For deformed nuclei with the doubly degenerate 
single-particle orbitals this approach supplemented by the appropriate use 
of symmetries and truncations was already shown to be quite effective |511 
E2]. The diagonalization of the general pairing Hamiltonian (|3J) is much 
simpler than that of the full shell-model Hamiltonian due to the possibility of 
classifying many-body states within the seniority scheme [3J HUSH], especially 
in the case of spherical symmetry (jSJ). Long ago it was shown that this 

approach is useful not only in the exactly solvable degenerate model but in a 
realistic shell model context as well. With a perspective to complement the 
pairing problem with the subsequent account for other parts of the residual 
interaction, we consider this path promising and quite practical. 

It is well known that the pair annihilation, Lj, pair creation Lj, and 
occupation number operator (shifted to the middle of the j-subshell), 



[Lj, L°,) = 6 jf L v [Lj, L° f ) = -S jf L], [Lj, L)] = -26j r L°. (26) 



Therefore the pairing Hamiltonian Q preserves all partial quasispins Aj, 



(25) 



form an SU(2) algebra of "quasispin" for each j-subshell, 



L 2 j = (L°) 2 + \{L)Lj + LjLj) = Aj(Aj + 1). 



(27) 



The partial seniority quantum numbers, 




2A 



(28) 
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are also conserved. They express the number of unpaired, and therefore not 
participating in the pairing interaction (J3|L particles. The fully paired j -level 
corresponds to the maximum partial quasispin Aj — fi/4 and lowest partial 
seniority Sj — 0. 

The pair transfer L'^Lj between the levels j — > j' changes the occupancies, 
i.e. projections L° and L°,, keeping intact the lengths of quasispins Aj and 
Aj/, and, whence, seniorities Sj and sy. The space is decomposed into sectors 
with given partial seniorities Sj, and the basis states within each sector can 
be labelled by the set of occupancies Nj under a constraint . Nj = N, 
the total valence particle number. The passive (unpaired) particles occupy 
fixed orbitals and create nonzero seniorities. They influence the dynamics 
indirectly, through the Pauli blocking. The states with zero total seniority 
s = Y2j s j have the total spin J = 0, while for s ^ the further decomposition 
with respect to the rotation group is possible, and some many-body states 
with different angular momentum coupling but the same seniorities remain 
degenerate. 

Using the states with given values of Sj and various possible occupancies 
Nj as a basis, it is easy to construct the Hamiltonian matrix that is essentially 
the matrix with respect to the sets of Nj. The diagonal matrix elements are 



Nj- C f(Nj-Sj)(n j -Sj-Nj + 2) 



({sj},{N j }\H p \{s j },{N j }) = Yl 

(29) 

Each term in the square brackets gives a full solution for the pairing problem 
on a degenerate j-level. Clearly, as long as seniority is small, Sj <C Qj, 
each unpaired particle increases energy by Aj = Gjjflj/4, and this quantity 
plays the role analogous to that of the energy gap in the BCS theory. The 
off-diagonal matrix elements for the pair transfer j' — > j are 

Hj^j = [ {Nf _ Si ,)(n i# _ Sjl _ at., + 2 )(Oj - 8i - N 3 )(N, - Sj + 2)] 1/2 

(30) 

The highest matrix dimension is encountered for the lowest possible total 
seniority, s = for an even number of particles, and s = 1 for an odd number 
of particles. But, even for heavy nuclei, this dimension does not exceed few 
thousand (in modern shell model computations one has to deal with dimen- 
sions 10 8 and higher in the m-scheme). In addition, the Hamiltonian matrix 
is very sparse. As a result of the numerical diagonalization we obtain the 
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d-5/2 


d-3/2 


Sl/2 


hn/2 




-6.121 


-5.508 


-3.749 


-3.891 


-3.778 


97/2 


0.9850 


0.5711 


0.5184 


0.2920 


1.1454 


^5/2 




0.7063 


0.9056 


0.3456 


0.9546 


d-3/2 






0.4063 


0.3515 


0.6102 


Sl/2 








0.7244 


0.4265 


hn/2 










1.0599 



Table 1: Single-particle energies and pairing matrix elements Vo(jj;j'j') (in 
MeV) for the shell model space from 100 Sn to 132 Sn; matrix elements are 
determined from G-matrix calculations. 

spectrum of states for a given set of seniorities. For example, for an even 
system, the condition s = selects all zero partial seniorities, Sj = 0. All 
those states correspond to pair condensates that differ by the distribution 
of the average partial occupancies {Nj} among the subshells. In a standard 
language of the BCS theory supplemented by the random phase approxima- 
tion (RPA), the excited states for s = are various pair vibrations. However, 
here we do not make any assumptions of boson character or harmonic spec- 
trum of excitations. The next section illustrates the typical results of the 
diagonalization. 

4 Example: A chain of even isotopes 

The longest known chain of tin isotopes is a subject of extensive experimental 
and theoretical studies. Even considering the proton subsystem, Z = 50, as 
an inert core, we have to deal with the neutron model space that is too 
large for a direct diagonalization. Modern computational techniques that 
use the Lanczos iteration method allow for exact determination of few low- 
lying states in systems with up to 12 valence particles (35] • These results are 
essential for testing the approximate techniques. It is known that pairing 
correlations play a major role in forming the ground state wave functions of 
tin isotopes. 

Unlike in many other nuclear systems, pairing in tin isotopes is quite 
strong and stable being sufficiently above the point of the BCS phase tran- 
sition. There is only a relatively minor weakening in the mass region near 
114 Sn due to a gap between d 5 / 2 and g 7 / 2 and the rest of single-particle or- 
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bitals. We specifically would like to explore this region in order to discuss 
physics of the BCS phase transition. 

For tin isotopes ranging in mass number from A = 100 to A = 132, we 
assume a configuration space between the two neutron magic numbers, 50 
and 82. The valence neutron space contains here five single-particle orbitals, 
hn/2, d-3/2, su2, gr/2, and d 5 / 2 ■ We adopt parameters shown in Table [TJ the 
single-particle energies taken from experimental data and the interaction ma- 
trix elements from the G- matrix calculation [33] . The interaction parameters 
Vb in Table 1 are related to the pairing strengths Gjj> as 

G S r = \vo(jj;j'j'W(2j + l)(2j' + l). (31) 

The shell model calculations with these parameters reproduce the spec- 
troscopy of tin isotopes in the region A = 120 to A = 130 with a good 
accuracy. In parallel we discuss similar effects in calcium isotopes where we 
used a well established FPD6 interaction [21] . The fp neutron valence space 
covers calcium isotopes from 40 Ca to 60 Ca. The weakening of pairing in Ca 
occurs near 48 Ca, i.e. for the f 7 / 2 subshell closure. Results of the calculation 
for the 114 Sn region are shown in Tableland for the region 48 Ca in Table [3] 

An important consequence of the proximity to the BCS phase transi- 
tion is a reduction of pairing correlation energy predicted by the BCS but 
not confirmed by the exact solution. In the tin example the BCS under- 
predicts binding energy by about 0.4 MeV, while for calcium this number 
reaches 0.6 MeV. A similar difference appears in one-nucleon separation en- 
ergy Sj(N) = E(N — 1) — E(N) (index j here denotes the orbital of the 
unpaired nucleon in an odd-iV system), as can be seen from the comparison 
of rows (e) and (f) with the corresponding BCS prediction, lines (m) and (n), 
Tables 2 and 3. As stressed in Ref. [21], this discrepancy can be crucial for 
the nuclei near drip lines. Unfortunately, the BCS hardly can be improved 
with respect to the treatment of weak pairing. Even for the complicated 
particle number projection techniques accompanied by the the variational 
procedures on a broader set of mean field states, it remains unclear to what 
extent it is possible to describe the pairing phase transition, and whether 
the high-lying pair vibrations are included, the step needed to account for 
missing correlation energy. 

Another related feature is the difference in predicted occupation numbers 
that can be inferred from comparing rows (a) and (b) with (i) and (j) in 
Tables 2 and 3. A proper account for this difference can partially help to 
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ff P ■ 

Hi r . 




(aj 
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6.96 


4.46 


0.627 


0.356 


1.6 


(b) 


n (N\ 
u](N) 
v](N) 


0.870 


0.744 


0.157 


0.178 


0.133 


(c) 


0.128 


0.252 


0.838 


0.817 


0.863 


(d) 


0.865 


0.736 


0.155 


0.177 


0.131 


(e) 


Sj(N+l) 


2.8 


3.13 


3.14 


3.39 


3.29 


(f) 


Sj(N) 


6.86 


6.55 


7.25 


6.98 


7.12 


(g) 


\(N + 2;0\P]\N; 0)| 


0.68 


0.779 


0.617 


0.514 


1.03 


00 


\{N;0\Pj\N-2; 0)| 


0.81 


0.93 


0.524 


0.396 


0.845 




fiCS : 




(i) 


TV, (TV) 


6.71 


4.14 


0.726 


0.507 


1.91 


(j) 


n,(iV) 


0.839 


0.69 


0.181 


0.254 


0.159 


(k) 


A, 


1.31 


1.43 


1.43 


1.38 


1.25 


(1) 


^•(ivV 1) 


1.78 


1.55 


1.86 


1.59 


1.71 


(m) 


2.89 


3.21 


3.11 


3.21 


3.26 


(n) 


Sj(N) 


6.89 


6.64 


7.2 


7.03 


7.06 


(o) 


bcs(0|P,|0)bcs 


0.734 


0.801 


0.545 


0.435 


0.896 



EP: £( 114 Sn) = -86.308, £( 116 Sn) = -95.942, E( n2 Sn) = -75.831. 
BCS: £( 114 Sn) = -85.938, fx = -5.035. 

Table 2: The results of the exact pairing solution (EP) compared to the 
BCS solution for the 114 Sn nucleus; the interaction matrix elements are de- 
termined by the G-matrix calculations, see text; the separation energies Sj, 
quasiparticle energies tj and the pairing gaps Aj are given in MeV. 
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BCS : 



BCS 



Nj(N) 
nj (N) 
A, 

Sj(N+ 1) 



BCS 



7/2 
-8.39 



3/2 
-6.5 



5/2 
-1.9 



6.87 0.85 0.173 

0.858 0.212 0.0288 

0.133 0.779 0.97 

0.848 0.212 0.0281 

4.48 5.78 1.55 

9.64 9.75 14.1 

0.706 0.928 0.289 

1.07 0.612 0.288 



1/2 
-4.48 

0.111 
0.0557 
0.939 
0.0555 
4.09 
11.5 
0.309 
0.232 



6.5 1.22 0.155 0.124 

0.813 0.304 0.0258 0.062 

1.66 1.44 1.73 1.53 

2.13 1.56 5.45 3.18 

4.5 5.66 1.64 4.1 

9.64 9.51 13.8 11.3 

0.78 0.651 0.275 0.241 



EP: £( 48 Ca) = -71.215, E( 50 Ca) = -85.149, £( 46 Ca) = -55.501. 
BCS: £( 48 Ca) = -70.591, /i = -7.335. 

Table 3: Properties of the weakly-paired 48 Ca nucleus; the FPD6 interaction 
was used in these calculations; all energies are given in MeV. 
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correct the binding energy. In the presence of additional interactions, the 
monopole contribution to the energy can be particularly sensitive to the 
precise occupation numbers. Furthermore, in the use of mean-field methods 
for paired systems, a good reconstruction of the density matrix generated by 
the pairing is of critical importance. 

The exact pairing treatment (EP) becomes increasingly important in con- 
sidering the reaction amplitudes with paired nuclei. The one-nucleon transi- 
tion amplitudes defined in the exact solution via Eq. (|13|) can be compared 
with the corresponding BCS quantities. Since, similarly to the recursive ap- 
proach, these amplitudes connect different nuclei, the standard BCS relations 
Vj = rij, u 2 j = 1 — rij and Uj+v? = 1 are no longer true. Deviations from these 
equalities are clearly enhanced in the phase transition region where adding 
an extra particle can make a sharp difference. The BCS theory with an un- 
certain particle number does not account for such effects. The pair emission 
amplitudes exhibit even larger differences. Rows (g) and (h) of the Tables 
2 and 3 show these amplitudes for adjacent even systems. The numbers are 
noticeably different (in the BCS approach they are substituted with a single 
set shown in row (o)). These discrepancies are particularly crucial for weakly 
bound nuclei since not only the binding energy is affected by the improved 
treatment of pairing but also there are significant corrections to the reaction 
amplitudes. 

Further insight into the situation can be gained by varying the coupling 
strength. For this purpose we consider a "ladder" model that contains ten 
double-degenerate single-particle orbitals equally spaced with the interval of a 
unit of energy. The valence space is assumed to be half-occupied with N = 10 
particles. The most interesting region is near the Fermi surface, eF = 5. For 
Fig. consider the first single-particle level above the Fermi energy. As 

in the previous example, the BCS result significantly deviates from the exact 
solution near the phase transition, around G = 0.5, as seen from panel (a). 
In the same region one can observe a slight difference between u|, 1 — 
and rij in the exact solution. For the pair emission process, the differences 
between the BCS and exact solution become more pronounced. Here the 
particle number uncertainty is crucial since the level under consideration is 
above Fermi energy for iV = 10, but below it for N = 12. 

Contrasting the exact solution with the mean field picture, we can notice 
that the occupation number operators Nj in general may have nonzero off- 
diagonal matrix elements between states of the same seniority. In Figs. [2J 
114 Sn, and|21 the ladder model, the matrix elements between the ground state 
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Figure 1: Comparison of BCS and EP solutions for the ladder system as 
a function of pairing strength G. In panel (a) occupation numbers for the 
sixth level (first level above the Fermi surface) and spectroscopic factors for 
capture and decay are compared to exact calculation. For the same level, 
the pair emission and pair absorption amplitudes are compared, panel (b), 
to the BCS prediction. 
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Figure 2: Matrix elements of the operator Nj for the level in 114 Sn 
between the ground state and all pair vibrational states (zero seniority) 
plotted as a function of excitation energy of the state 0. 



18 



and all s = states are shown as a function of excitation energy. In all cases 
the off-diagonal matrix elements rapidly fall off. In the case of weak pairing 
(Fig. El circles), one can still see the structure of excited states based on 
the equidistant single-particle spectrum. For stronger pairing (compared to 
the single-particle level spacing), Fig. EH triangles, as well as in the realistic 
case for spherical symmetry, Fig. |21 the decrease of matrix elements is more 
uniform and can be approximated in average by an exponential function of 
excitation energy. This indicates chaotization of motion even in the sector 
with seniority s — 0, [SJ]]. Only very few states with relatively large matrix 
elements may carry pair-vibrational features. As follows from the extended 
shell model analysis |57|, the exponential tails of the strength functions are 
typical for many-body quantum chaos [HE] • The property of exponential con- 
vergence was demonstrated in Ref. j2H] and the extrapolation based on this 
property was later used [23 EM EH] as a practical tool for getting reliable 
quantitative results in shell-model calculations of intractable large dimen- 
sions. 



5 Continuum effects 

As the main interest of low-energy nuclear physics is moving to nuclei far 
from the valley of stability, the continuum effects become exceedingly im- 
portant for a unified description of the structure of barely stable nuclei and 
corresponding nuclear reactions. The pairing part of the residual interaction 
in some cases is the main source of the nuclear binding; spectroscopic factors 
and reaction amplitudes are also critically dependent on pairing. 

As a demonstration of a realistic shell-model calculation combining the 
discrete spectrum and the continuum, we consider oxygen isotopes in the 
mass region A = 16 to 28. In this study we use a universal sci-shell model 
description with the semi-empirical effective interaction (USD) [23]. The 
model space includes three single-particle orbitals lsi/2, 0d 5 / 2 and 0d 3 / 2 with 
corresponding single-particle energies —3.16, —3.95 and 1.65 MeV. The resid- 
ual interaction is defined in the most general form with the aid of a set of 63 
reduced two-body matrix elements in pair channels with angular momentum 
L and isospin T, ((J3T3, j&T4)LT\V\{jiT\, j2T 2 )LT), that scale with nuclear 
mass as (A/18) -0 - 3 . 

In the discrete spectrum the full shell model treatment is possible for 
such light systems. Aiming at the study of the continuum effects, that in- 
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Figure 3: Matrix elements of N 6 (the sixth single-particle state in the ladder 
model) between the ground state and excited s = states as a function of 
excitation energy. 



20 



-8 MeV 




-24 MeV 



-32 MeV 



-40 MeV 




Figure 4: Schematic picture of lowest seniority states and possible decays 
in oxygen isotopes. Dotted arrows indicate the decays blocked by seniority 
conservation. 

crease significantly the computational load, here we truncate the shell-model 
space to include only seniority s = and s — 1 states. This method, "exact 
pairing + monopole", is known j2H] to work well for shell model systems 
involving only one type of nucleons (in the case of the oxygen isotope chain 
only neutrons and the interaction matrix elements with isospin T = 1 are in- 
volved). The two important ingredients of residual nuclear forces are treated 
by this method exactly: the monopole interaction that governs the binding 
energy behavior throughout the mass region, and pairing that is responsible 
for the emergence of the pair condensate, renormalization of single-particle 
properties and collective pair vibrations. 

In the resulting simplified shell-model description, the set of the original 
30 two-body matrix elements in the T = 1 channel is reduced to 12 most 
important linear combinations. Six of these are the two-body matrix elements 
for pair scattering in the L = channel describing pairing, and the other six 
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correspond to the monopole force in the particle-hole channel, 

V j]r = Y,{lL+l)Uj')Ll\VUj')Ll), (32) 

where j and j' refer to one of the three single-particle levels. 

We assume here that the 0<i3/2 orbital belongs to the single-particle con- 
tinuum and therefore its energy has an imaginary part. In this model we 
account for two possible decay channels |c) for each initial state |$), a one- 
body channel, c = 1, and a two-body channel, c = 2. The one-body decay 
changes the seniority of the 0d 3 /2 orbital from 1 to in the decay of an 
odd- A nucleus and from to 1 for an even- A nucleus. The two-body decay 
with the zero angular momentum of the pair removes two paired particles 
and does not change the seniority. The two channels lead to the lowest en- 
ergy state of allowed seniority in the daughter nucleus; transitions to excited 
pair-vibrational states are ignored. This results in 

es/a(S) = e 3/2 - ^3/2 (E* - £ (1) ) 5/2 - 1 a 3/2 (E* - E^f' 2 , (33) 

(c) 

where we assumed that one- and two-body decay parameters 7] are related 
as 7gy 2 = 73/2/2 = 73/2, and the particles are emitted in the <i-wave with 
£ = 2 . The energy dependence of the widths near decay thresholds E^ c ' is 
very important; 03/2 is the reduced width parameter that differs from 73/2 
by the absence of the energy factor. Three states with the valence particle 
at one of the single-particle orbitals can be identified as the 5/2 + ground 
state and l/2 + and 3/2 + excited states in the spectrum of 17 0. Their en- 
ergies relative to 16 correspond to the single-particle energies in the USD 
model. Experimental evidence indicates that the 3/2 + state decays via neu- 
tron emission with the width T( 17 0) = 96 keV. This information allows us 
to fix our parameter a 3/2 = r( 17 0)/(e 3/2 ) 5/2 = 0.028 (MeV)" 3 / 2 . Other two 
states are particle-bound, 71/2 = 75/2 = . 

With complex single-particle energies, the non-Hermitian effective Hamil- 
tonian for the many-body system is constructed in a regular way jHSES]- The 
Hamiltonian includes the Hermitian pairing and monopole terms as well as 
the energy-dependent non-Hermitian effective interaction through the open 
decay channels, the structure of which is dictated by unitarity jHSl- Energy 
dependence is determined by the proximity of thresholds, Eq. (J3~3*|) . We 
move along the chain of isotopes starting from 16 O. In this way, for each A 
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the properties of the possible daughter systems A — 1 and A — 2 are known. 
The chain of isotopes under consideration is shown in Fig. which includes 
s = and s — 1 states and indicates possible decays. The decays indicated 
by the dotted arrows are blocked in our model due to the exact seniority 
conservation. Non-pairing interactions in the full shell model mix seniori- 
ties making these decays possible. Since the effective Hamiltonian depends 
on energy and all threshold energies have to be determined self-consistently, 
we solve this extremely non-linear problem iteratively. We start from the 
shell-model energies -E s . m . corresponding to a non-decaying system. Then 
the diagonalization of the Hamiltonian at this energy allows us to deter- 
mine the next approximation to the complex energies £ = E — {i/2)T that 
give the position and the width of a resonance. This cycle is repeated until 
convergence that is usually achieved in less than ten iterations. 

The results of the calculations and comparison with experimental data for 
oxygen isotopes are shown in Table 0] Despite numerous oversimplifications 
related to seniority truncation, limitations on the configuration mixing and 
restrictions on possible decay channels and final states, the overall agreement 
observed in Tabled is quite good. If experimental data are not available, the 
results can be considered as predictions. The main merit of this calculation 
is in demonstrating the power and practicality of the EP method extended 
to continuum problems. The same calculation predicts also [HI] the cross 
sections of the processes related to the included channels providing the unified 
description of the structure and reactions with loosely bound nuclei. 

6 Thermal properties 
6.1 BCS approach 

The properties of the dense spectrum of highly excited states are usually 
described in statistical terms of level density, entropy and temperature. The 
shell model analysis ^§1 EH] revealed a certain similarity between many- 
body quantum chaos and thermalization. In particular, the Fermi-liquid 
approach to the complex many-body system modelled as a gas of interacting 
quasiparticles, turns out to be applicable not only in the vicinity of the Fermi 
surface but even at high excitation energy. 

Here we consider the thermalization properties of the paired system. The 
BCS operates with the quasiparticle thermal ensemble. The expectation 
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A 


J 


E (MeV) 


T (keV) 


E cxp . (MeV) 


r cxp . (keV) 


16 





0.00 





0.00 





17 


5/2 


-3.94 





-4.14 





17 


1/2 


0.78 





0.87 





17 


3/2 


5.59 


96 


5.08 


96 


18 





-12.17 





-12.19 
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5/2 


-15.75 





-16.14 





19 


1/2 


1.33 





1.47 





19 


3/2 


5.22 


101 


6.12 


110 
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-23.41 





-23.75 
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5/2 


-26.67 





-27.55 
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1/2 


1.38 
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3/2 


4.60 


63 
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-33.94 





-34.40 
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1/2 


-35.78 





-37.15 
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2.12 
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2.57 
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-40.54 





-40.85 





25 


3/2 


-39.82 


14 






25 


1/2 


2.37 









25 


5/2 


4.98 
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-42.04 









27 


3/2 


-40.29 


339 
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1/2 


3.42 


59 






27 


5/2 


6.45 


223 






28 





-41.26 


121 







Table 4: Seniority s = and s = 1 states in oxygen isotopes. Energies and 
neutron decay widths are shown. Results are compared to the known data. 
Ground state energies relative to the 16 O core are given in bold. The rest of 
the energies are excitation energies in a given nucleus. 
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value for an occupancy of a given single-particle state is 

n = (a^a) = u 2 {ara) + v 2 {aa { ) + uv(arafi + oia). (34) 

Under the assumption of the thermal equilibrium quasiparticle distribution, 
the last term in ()34j) disappears, while the first and the second terms give 

(ojaj) = uj, (dcjOj) = l — v~-. (35) 

The occupation numbers for quasiparticles are defined by the Fermi distri- 
bution with the zero chemical potential, 

u J (T) = [l + exp(Pe j )}- 1 , (3 = 1/T, (36) 

and temperature-dependent quasiparticle energies 



e J (T) = ^(e J -^ + A](T), (37) 

The thermal evolution of the pairing gap A(T) is determined by the self- 
consistent BCS equation with the quasiparticle blocking factor included, 

^ [1 - Vf{T) — vr,(T)]Af(T) 
^ = E G " 2e f (T) ■ ^ 

The occupation numbers of original particles are given by 

rij = u 2 Vj + v 2 (l — Uj), (39) 

where the coherence factors u and v also depend on temperature via e 3 -(T). 
(This discussion is closely related to Ref. jOSl-) As a result, we obtain 
equations for the gap (|55|) and chemical potential using Eq. ()39p at a given 
external temperature T that governs the quasiparticle distribution. 



6.2 Statistical spectroscopy of pairing 

A form of the pairing Hamiltonian allows for a relatively simple calculation 
of its spectroscopic moments. In this section we limit our consideration to 
the zero seniority block of a ladder system of total capacity Q with double- 
degenerate orbitals; the generalization for more realistic cases is straightfor- 
ward. In the ladder system the diagonal matrix elements simply renormalize 
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single-particle energies. It is convenient to set the chemical potential to zero 
and use variables £\ = e\ — Gn/2 following Eq. (|2L)p . We also denote the 
off-diagonal pairing matrix elements as Q12 = (1 — 5\2)Gi2- 

The centroid of the distribution is determined by the single-particle spec- 
trum, 

((E)) = Ne, (40) 

Here the double brackets imply averaging over all many-body states, while 
the overline means averaging over single-particle states according to the def- 
inition 

G* = ^Tr(g k ). (41) 

The second moment of the distribution, the variance, is a sum in quadratures 
of the single-particle width and the width due to pairing, 

a 2 = (((E - ((E))) 2 )) = 2N ^J 2 N) • (42) 

The third moment, the skewness, indicates deviations from the normal dis- 
tribution. It is given by 

All odd central moments are asymmetric in Q and thus vanish for Q = 0. 
The skewness is also a special case since in the ladder system it does not 
depend on single-particle energies. For attractive pairing the skewness is 
always negative indicating a longer tail of the distribution towards the lower 
energies. This supports the pairing character of the low-lying states that due 
to their collective nature are pushed further down from the centroid of the 
distribution. 

The density of states p{E) allows for a thermodynamical determination 
of the temperature, 

f = ^(p(E)). (44) 

Despite the presence of higher moments, the density of states of paired sys- 
tems, as in a more general class of two-body Hamiltonians [T^j, can be closely 
approximated by the Gaussian distribution. An actual distribution p(E) is 
shown in Fig. Assuming the Gaussian distribution with the mean value 
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(f4"U|) and the width a, Eq. (fi"2J). the temperature as a function of energy can 
be found [IH] using flUj) , 



a 2 



T ™ = J^e- (* ) 

The negative T branch is an artifact of the finite Hilbert space. 

The Gaussian distribution is getting distorted when the pairing becomes 
strong and the low-lying states very collective. A minor manifestation of this 
collectivity is seen in Fig. 0(b) for G = 1. As G grows, the deviations from 
the Gaussian shape become more transparent clearly revealing the seniority 
structure as seen in Fig. |H1 



6.3 Quasiparticle temperature 

As was discussed in detail in |19| , there is no unique definition of tem- 
perature in a self-sustaining isolated mesoscopic system. Complementary to 
the thermodynamic (microcanonical) definition of the previous subsection, 
we can find the effective value of temperature for each individual many-body 
state by fitting the occupation numbers found in the EP solution to those 
given by the thermal ensemble, Eq. We can refer to this assignment 

as a measurement with the aid of the quasiparticle thermometer, and denote 
the resulting temperature as T . 

The correspondence between excitation energy and quasiparticle temper- 
ature T for each eigenstate in the 12-level model is presented in Fig. The 
scattered points clearly display a regular trend to thermalization in agree- 
ment with the hyperbolas predicted by Eq. ([45)1 . Although thermodynamic 
and quasiparticle temperature are well correlated, Fig. |H1 their numerical 
scales are different, T « 2.5T. Furthermore, the concept of quasiparticle 
thermalization is meaningful only for relatively weak pairing, where large 
fluctuations due to the proximity of the phase transition are present. The 
quality of thermalization deteriorates as stronger pairing makes the dynamics 
more and more regular. The inability of the strongly paired system to fully 
thermalize dynamics was demonstrated earlier (HE]. The role of non-pairing 
interactions is essential for equilibration. But the failure of the single-particle 
thermometer to reflect correctly the spectral evolution in the limit of very 
strong interaction is a general feature [T§1 lol)] . 
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Figure 5: Density of states as a function of energy is shown by the histogram. 
The solid curve indicates the Gaussian distribution with the parameters de- 
termined by Eqs. (J4"U|) and ([42)1 . Both panels refer to half-occupied systems. 
Panel (a) for the 12- level system with G = 0.5 displays an overall good Gaus- 
sian fit while the deviations are more clear in panel (b) where the density of 
states is plotted in the log scale for a larger system with stronger pairing, 
G — 1. Here the effect of negative skewness is visible, as well as the presence 
of the gap. 



28 




Figure 6: Density of states in the half-occupied system with f2 = 28 and pair- 
ing strength G — 3. The strong pairing dominates single-particle structure 
resulting in the peaks of the level density that correspond to the quasispin 
families of the degenerate model. 
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Figure 7: Excitation energy and quasiparticle temperature T of individual 
states in the half-occupied, 12-level ladder model. The solid line shows the 
energy-temperature relation from Eq. (J45|) . assuming scaling T « 2.5T. 
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Figure 8: Comparison of thermodynamic and quasiparticle temperatures for 
the half-occupied, 12-level ladder model. Three panels correspond to pairing 
strengths G =0.5, 0.7 and 1.0, from left to right, respectively. The regular 
trends in scattered points confirm approximate thermalization with temper- 
ature ratio T « 2.5T. The agreement deteriorates for stronger pairing. 
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6.4 Pairing phase transition 

In Fig. Elthe pairing gap is plotted as a function of quasiparticle temperature. 
The gap was calculated using Eq. (38) with the quasiparticle temperature 
substituting the BCS external temperature parameter. This makes the con- 
sideration consistent with the occupancies given by Eqs. (36) and (39). The 
half-occupied, 12-level ladder model was again used for this example. The 
choice of a larger system not only results in the increased number of s = 
states, but, more importantly, reduces the particle number fluctuations that 
can disrupt the fitting procedure, especially in the pairing phase transition 
region. 

Fig. demonstrates the phase transition from the paired state at lower 
temperature (or excitation energy) to a normal state at a higher temperature. 
Few low-lying states have a considerable pairing gap, whereas the gap disap- 
pears in the sufficiently excited states. The invariant entropy jHZI can be an 
alternative method for visualizing the phase transition 168) . This quantity 
is basis-independent and reflects the sensitivity of a particular eigenstate to 
the changes of a parameter of the many-body Hamiltonian, here the pairing 
strength G. The peak in the invariant entropy points out the location of 
the pairing phase transition as a function of G and excitation energy E of a 
particular state. However, it is important to notice that the phase transition 
pattern is strongly influenced ^21 EH] by other parts of residual interaction. 

7 New theoretical perspectives 

7.1 Hartree-Fock approximation based on the exact 
pairing solution 

Instead of the normal Fermi-occupation picture with a Slater determinant as 
a trial many-body function for the mean field approximation, the EP solution 
with its specific single-particle occupancies provides a new starting point for 
the consistent consideration of other parts of the residual interaction. In this 
subsection we illustrate this point with the help of Belyaev's |6R 4 pairing plus 
quadrupole (P+Q) model Hamiltonian for a single- j level, 




(46) 
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Figure 9: Pairing gap as a function of quasiparticle temperature in the 12- 
level model with pairing strength 1^=0.5 and 0.7 for left and right panels, 
respectively. 
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where the multipole operators are defined as 

m\m2 

Only the ratio xjG of the strength of quadrupole-quadrupole interaction to 
the pairing strength is important since the energy scale can be fixed so that 
G = 1. 

In the pure pairing limit, \ = 0, the degenerate pairing model is recovered 
with the ground state energy (Q = 2j + 1) 

E = G^-(tt- N + 2). (48) 

Pairing correlation energy in the BCS with constant pairing is given by 
A 2 /G . For the exact solution we define the correlation energy as the ground 
state expectation value of the pairing part of the Hamiltonian with the 
monopole contribution subtracted. In the degenerate case 

N 

E COII = G-(n-N). (49) 

The opposite limit with no pairing can be treated by making a transition 
to a deformed mean field in the Hartree approximation [70]. For axially 
symmetric deformation, the expectation value of the quadrupole moment is 

. . 2[3m 2 - j(j + 1)1 , , 

(M 20 =V 7= ===^= n m , 50 

in terms of the occupation numbers n m = (o^ a m ) in the intrinsic frame with 
the z-axis oriented along the symmetry axis. In this case 

(M2-2) = (M22) = 0. (51) 

The deformed single-particle energies in the body-fixed frame can be obtained 
via the usual self-consistency requirement, 

em -~ X ^-l)(W-2) (M20) - (52) 

The energy minimum for an even-iV system corresponds to the Fermi occu- 
pation of the N/2 lowest pairwise degenerate orbitals \m\ = 1/2, 3/2, . . . (iV — 
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l)/2 for prolate, or \m\ = j, j — 1, . . . j — (N — 2)/2 for oblate shape. The 
corresponding quadrupole moment is then given by 

111 \ 1 N(Q 2 -N 2 ) 

(M 2 o) = ; (53) 

for prolate deformation and 

1 N(2Q - JV)(ii - N) 
{Mm) Vii(ii^l)(ff-4) (M) 

for oblate deformation. Deformation energy is quadratic in so that 

the oblate deformation is preferred for N < Q/2 and the prolate one for 
N > Q/2. 

The full problem is driven by the competition between pairing and defor- 
mation. While deformation tends to split the single particle energies (Nilsson 
orbitals), Eq. (|52jl . the pairing can resist such a shape transition by creating 
a particle distribution unfavorable for deformation. In the single- j model 
these effects have been discussed by Baranger [70] with the help of the BCS 
and the Hartree approximation. However, as demonstrated earlier, the BCS 
may be unreliable in the transitional region. The use of the projected HFB 
and the resulting improvement against traditional HFB for a similar single- j 
model have been recently discussed in [7_T]. Our goal here is to supplement 
the Hartree treatment with the exact pairing solution. 

Similar to the Hartree+BCS approach in [70] we look for a self-consistent 
solution, where the occupation numbers in the deformed basis agree with 
the exact solution to the pairing problem (22] . Prior to calculations, one can 
estimate the ratio x/ G corresponding to the BCS phase transition. Assuming 
for example an oblate shape with iV < Q/2 we have the Fermi energy at 
m = j — (N — 2)/2 with the density of single-particle states approximately 
found as u F = 2/(e m _i — e m ), that for a half-occupied system leads to 

(X/G)«it»3.6n. (55) 

In Fig. El we present the results for a model with j = 19/2, The particle 
number, N = 8, was selected to avoid an exactly half-occupied shell when 
particle-hole symmetry and oblate to prolate shape change lead to special 
features that are of no interest for our goal. For this model the transition from 
spherical to deformed shape takes place at x/G ~ 70, Eq. (|55[). The pairing 
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Figure 10: Comparison of pairing correlation energy, panel (a), and 
quadrupole deformation, panel (b), for Hartree+BCS and Hartree+EP treat- 
ments of the P+Q model for a single- j level, j = 19/2, and N = 8 particles. 
For pairing correlation energy G = 1.0 is assumed. 



correlation energy shown in Fig. ITUT a) as a function of xjG starts near 
x/G = with a value prescribed by the degenerate model. As the relative 
strength of the quadrupole-quadrupole interaction increases, the deformation 
inhibits pairing. However, in Fig. [TUT a) we see a key difference between the 
BCS and exact treatment. Within the BCS the pairing correlation energy 
goes to zero quite sharply once the system becomes deformed. In contrast, 
the exact solution finds that pairing correlations decay very slowly and extend 
far into the deformed region. In panel ITW b) the expectation value of (Ai 2 o) is 
shown as a function of xjG. In the pairing limit x/G — > the deformation is 
zero whereas in the deformed limit the value (.M20) is expressed via Eq. ()54|) 
or (jSHJ). Here again the exact treatment produces a softened and extended 
phase transition. 
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Figure 11: Comparison of pairing correlation energy and quadrupole de- 
formation; in contrast to Fig. El the exchange contribution from the 
quadrupole-quadrupole interaction to the pairing channel is included. 

The Hartree approximation ignores another important effect relevant to 
our consideration, namely, the contribution of the exchange terms to the pair- 
ing channel. This contribution is particularly strong in small systems and 
leads to an additional enhancement of the pairing strength G —>■ G + 2x/f2 2 
[73] . The results of Hartree- Fock + EP calculations that include this addi- 
tional term are shown in Fig. ^2 E)ue to the exchange term, pairing correla- 
tions never disappear. Presence of pairing correlations in a pure quadrupole- 
quadrupole Hamiltonian is also confirmed by an exact solution in the full shell 
model diagonalization |73| . The BCS treatment, however, fails to reproduce 
this effect. 
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7.2 New Random Phase Approximation 

Here we show how the RPA-like approximation can be developed starting 
from the exact solution of the pairing problem. There are two main types 
of the RPA used in the literature (and a variety of close approaches distin- 
guished by the details of the formalism), the RPA based on the vacuum of 
noninteracting particles or that of the BCS, or HFB, quasiparticles (the so- 
called QRPA). We will try to describe collective vibrations generated by the 
residual interaction on top of the exact ground state of the pairing problem. 
The formalism of the generalized density matrix [TJJ [75] seems to be suitable 
for this purpose. 

The generalized density matrix (GDM) is the set of the operators 

R\2 = d\, a i (56) 

acting in the full Hilbert space of a many-body system; this set at the 
same time forms a matrix labeled by single-particle subscripts {1,2). We 
do not perform any canonical transformation to the quasiparticle operators 
and therefore work invariably within a system of a certain particle number. 
The one-body observables as operators in many-body space are traces of the 
GDM operator over single-particle indices, 

Q = ^2% Q = ^2 ?ia a i a 2 = Tr (Q R )- ( 57 ) 

a 12 

With the Hamiltonian of the system taken as a sum of independent particle 
energies in the mean field, e±, and the general residual two-body interaction 
^i2;34j the exact operator equations of motion for the GDM can be symboli- 
cally written as 

[R,H] = [S,R], (58) 

where S is the generalized self-consistent field operator (a linear functional 
of the GDM), 

S = e + W{R}, W u {R} = J2Vi2;3±R 3 2, (59) 

23 

and the interaction matrix elements are antisymmetrized. 

Now we assume that the Hamiltonian contains the pairing part (J1J) as 
well as other residual interactions, H = H p + H' . Correspondingly we can 
put 

R = R° + R', W = W° + W, W = W{R°}, W' = W{R'}. (60) 
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The assumption of the exact solution of the pairing problem means that we 
found the occupancies, R°, and the pairing potential, W°, satisfying 

[R ,H p } = [e + W°,R°}. (61) 

This stage of the solution provides the states \s, a) with seniority s and energy 
E sa , where a numbers the states within the subset of certain seniority; if 
needed we also can explicitly indicate rotational quantum numbers J a and 
M a . The remaining part R' of the GDM should satisfy 

[R', H p ] + [R°, H'] + [R', H'] = [W, R°] + [e + W°, R') + [W, R'}. (62) 

The commutators in such expressions are to be understood as, for example 

[W, R'] 12 = J2( W ™ R v - ^3^32)- (63) 

3 

This is the point where we can make RPA-like approximations. 

For definiteness, we consider the transitions from the paired states s = 
0, J = to the states with s = 2, J ^ in the next sector. We are looking 
for a collective mode that is related to such excitations. This means that 
there exist states, in our case coherent combinations of excited states with 
certain J, that have large off-diagonal matrix elements of excitation by a one- 
body multipole operator from the ground state. The latter can be in turn 
renormalized by the collective mode. Let us characterize this branch of the 
spectrum with the help of collective coordinates a and conjugate momenta 
7i (we omit in this symbolic derivation their quantum numbers of angular 
momentum and its projection). These variables are Hermitian quantum op- 
erators that satisfy the commutation relation [a, it] = i so that no procedure 
of subsequent requantization is needed. The collective Hamiltonian of the 
mode can be written as 

H' = \cc? + ^+ (64) 

where the scalar contraction of the tensor operators is implied and the dots 
include high order anharmonic terms important for the soft mode [7d1 I77j . 

We are looking for the operator solution of Eq. (162(1 in the form of the 
expansion in collective operators a and ir, 

R' 12 = r^a + r^n+ (65) 
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and 

W{ 2 =wg 0) a + wf ) 2 1) n+ (66) 

where the superscripts (n, m) refer to the component containing n collective 
coordinate and m collective momenta operators. The dots again denote the 
higher order parts, n + m > 1, symmetrized in due way [76| 177] . The col- 
lective operators producing the transition in many-body space are written 
explicitly whereas the coefficients rj 2 m ' ) and w^, are the c-numbers to be 
found as single-particle amplitudes of the coherent superposition that forms 
a collective mode. The operator expansion (J551 IBljj) does not assume the 
smallness of anharmonic effects - we merely decompose the problem in var- 
ious operator structures. In the present context we limit ourselves by the 
harmonic part although it can be used EZj for the situations of strong 
anharmonicity as well. 

The operator R, by definition ([56)1 . has seniority selection rules |As| = 2. 
We take in the operator equation ()62j) matrix elements (0\c) between the 
ground state |0) and the collective state |c) in the adjacent sector s = 2 with 
angular momentum corresponding to that of collective operators a and ir. 
Now we evaluate the matrix elements of various terms in Eq. (J62)) aiming at 
the segregation of terms linear in a and 7r. The first term in the right hand 
side gives, according to Eq. (|B3|1. 

[R'i2, H p ] = (E° c - E°)(r^a + rffV). (67) 

Here the barred energies are the centroids of the energy distribution of the 
actual ground state and the one-phonon state in the sectors s = and s — 2, 
respectively. The second term in eq. (J62|) does not have required matrix ele- 
ments whereas in the commutator [R', H'] we need to perform commutation 
explicitly using the assumed form of the collective operators (|65H67|) . 

[m a ,H f \ = ^r^-iCrS ) a. (68) 

The situation with the terms and [W, R°] is more complicated. 

Within each sector of given s, the pairing solution GDM R° has not only 
diagonal but also off-diagonal elements between the eigenstates of H p . As 
shown in Figs. 2 and 3, very few pair- vibrational states have significant off- 
diagonal matrix elements of this type. For our illustrative purposes, here we 
neglect the off-diagonal terms within a given sector and take into account only 
diagonal elements of R° and W°. The neglected contributions correspond to 



40 



the anharmonic admixtures of pair vibrations to multipole modes and can 
be easily included in the consideration. Because of the specific character 
of the monopole pairing interaction, the matrix elements of R° and W° are 
diagonal over single-particle subscripts as well. Higher order structures in 
the collective Hamiltonian and in the GDM, as well as terms generated by 
the commutator [W, R'] do not contribute to matrix elements linear in a 
and 7r and with the selection rule As = 2. But, similarly to Eq. (J67j) . the 
commutators with R° and W° bring in the differences of the single-particle 
occupancies and pairing potentials averaged over the states contributing to 
the collective mode. 

As a result, we come to the coupled equations for the coordinate and 
momentum RPA amplitudes (these contributions can be distinguished by 
their behavior under time- reversal operation in the sector with J ^ 0): 

r^Q 12 - iCrU = [n 2 (c) - %(())]<, (69) 



rgOia + ^r{° 2 = [n 2 (c) - f»i(0)]«;g. (70) 

Here the generalized frequencies are introduced, 

n 12 = E° - E° + e 2 - Cl + W° c - W°. (71) 

This set of equations leads to a formal solution, analogous to that in the 
conventional RPA, 

io _ [n 2 (c) - wi(0 )] (n _ 10 , .^_ i01> 
- cu 2 



ri 2 ° = ^2 r J (Qu< + zC<), (72) 



-12 

C) 



ni _ Mc) - rii(O)] 01 % 10 



(fii2< + ^0, (73) 



-12 

where the unknown collective frequency is uj = (C/B) 1 ^ 2 . 

The collective elements of the GDM are to be found from the integral 
equations (fT2~j) and (|73)) with the specific choice of the residual interaction 
self-consistently generating the field W, Eq. (JBDJ). This can be done explicitly 
in the case of the factorizable multipole-multipole force; the frequency u and 
correct normalization of the mode are also obtained in this process similarly 
to the standard procedure in terms of the barred quantities. Having at 
our disposal the phonon amplitudes and we can self-consistently 
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find the barred quantities averaged over the collective wave functions. This 
procedure, that reminds the thermal RPA built on the equilibrium density 
matrix but with occupation numbers and mean-field corrections defined by 
the interaction rather than by an external heat bath, can be performed in an 
iterative manner. The cranking description for the deformed nucleus can be 
also reformulated in the same spirit; it is interesting to note that the pure 
EP solution predicts jHE] & n yrast-line with the moment of inertia close to 
the rigid body value. 

8 Conclusion 

The pioneering work by Belyaev jS] carried out a detailed analysis of pairing 
phenomena in nuclei. Applying the BCS techniques in the nuclear shell 
model environment, he demonstrated the effects of pairing on various nuclear 
properties, including the ground state structure, single-particle transitions, 
collective vibrations, onset of deformations and rotations. It was shown that 
the Cooper phenomenon in systems with a discrete single-particle spectrum 
does require, in contrast to large systems, a certain strength of the pairing 
interaction. The drawbacks of the BCS approximation, related to the particle 
number violation and the sharp disappearance of pairing correlations at the 
phase transition point, were also pointed out. 

The development started with Belyaev's work and supported by similar 
studies |7H1 EHl was continued through the next forty years. Now the pairing 
problem is well and alive being one of the main chapters of modern nuclear 
physics and mesoscopic physics in general. The interest in pairing is con- 
stantly revived by the accumulation of data and especially by the advances 
towards nuclei from stability where the pairing is a key tool that determines 
the binding of a system and its response to the excitation. At this point it 
becomes increasingly important to get rid of the shortcomings of the BCS 
approximation and unify the description of the structure and reactions. 

We presented a way of solving the pairing problem essentially along the 
lines similar to that of Belyaev's paper, substituting the BCS approximation 
by the exact solution simplified by the seniority symmetry. As a magnifying 
glass, this solution reveals and fixes the weak points of the standard approach. 
We saw the importance of the exact treatment for the ground state, low- 
lying excitations, coupling through the continuum, and spectroscopic factors 
associated with single-particle removal and pair emission (transfer) reactions. 
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We could also discuss on the new basis the global properties of the spectrum, 
thermalization and the phase transition region. In many cases, this exact 
treatment of pairing is practically simpler than solving the BCS equations 
with necessary corrections. 

Certainly, the pairing problem is only a part of the physics of strongly 
interacting self-sustaining systems. Other interactions, with their own coher- 
ent and chaotic features, should be included in the consideration. We gave 
preliminary answers to the questions of further approximations necessary for 
the cases when the full problem does not allow for a complete solution. New 
generalizations of the mean field approach and random phase approximations 
can be developed on the background of the exactly found paired state. The 
interplay of pairing and other residual interactions can be an exciting and 
practically important topic of future studies. 

The inspiring influence of Belyaev's ideas is gratefully acknowledged by the authors 
who belong to the first and third generation of his pupils. The work was supported 
by the NSF grant PHY-0070911 and by the U. S. Department of Energy, Nuclear 
Physics Division, under contract No. W-31-109-ENG-38. 
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